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^ . Abstract 

>. 
0^ . Violation of the Equivalence Principle (VEP) can lead to neutrino oscillation through the non- 

1_1/ ' diagonal coupling of neutrino flavor eigenstates with the gravitational field. The neutrino energy 

^\^ '■ dependence of this oscillation probability is different from that of the usual mass-mixing neutrino 

^P . oscillations. In this work we explore, in detail, the viability of the VEP hypothesis as a solution to 

the solar neutrino problem in a two generation scenario with both the active and sterile neutrino 

alternatives, choosing these states to be massless. To obtain the best-fit values of the oscillation 

parameters we perform a x^ analysis for the total rates of solar neutrinos seen at the Chlorine 

(Homestake), Gallium (Gallex and SAGE), Kamiokande, and SuperKamiokande (SK) experiments. 



o 
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Oh: 

n\ ' We find that the goodness of these fits is never satisfactory. It markedly improves, especially for VEP 

1-^ ■ transformation to sterile neutrinos, if the Chlorine result is excluded from the analysis. The 1117-day 

J> , SK data for recoil electron spectrum are also examined for signals of VEP oscillations. For these fits, 

we consider variations of the Standard Solar Model by allowing the absolute normalizations of the 



X 



^ ' ''B and hep neutrinos to vary. Here the fits are quite good but the best fit values of the parameters 
. . 1 are rather different from those from the total rates fits. A combined fit to the total rates and recoil 
electron spectrum data is also performed. We present the 90% confidence limit contours for all the 
three analyses mentioned above. The best-fit parameters obtained from the recoil electron spectrum 
and the combined analysis of rate and spectrum are used to predict the charge current and scattering 
electron spectrum at SNO. 
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I Introduction 

Oscillation of neutrinos from one flavor to another is currently the favored solution to the 
solar neutrino problem [0, |^, ^. This proposition has been strengthened by the recent 
SuperKamiokande (SK) atmospheric neutrino data Q which also support the existence of 
neutrino oscillations. The usual formulation of neutrino oscillations rests on two essential 
properties; namely, (a) the neutrinos are massive and are further not mass degenerate, and 
(b) the flavor eigenstates {i.e., u^-, ^^i-, and Vt) are not themselves the eigenstates of mass but 
rather linear superpositions of the latter. Distinct from the simplest vacuum oscillations, 
neutrino flavor conversion can also be induced by the passage of neutrinos in matter (in this 
case solar matter) due to the difference in the strength of weak interactions of neutrinos of 
different flavor with ambient electrons. This is the Mikheyev-Smirnov-Wolfenstein (MSW) 
effect. Vacuum oscillations and the MSW effect are widely consider d to be strong 
candidates for the solution of the solar neutrino problem. Much work has already been 
done to analyze the available solar neutrino data in terms of these alternate possibilities 
and substantial effort is still being devoted to obtain the parameters that best fit the data 
§0,1,1,0,0. 

Though the vacuum oscillation and MSW solutions relying on massive neutrinos (we 
refer to these collectively as 'mass-mixing' solutions henceforth) are by far the most popu- 
lar scenarios for addressing the solar and atmospheric neutrino data, oscillations can also 
originate from other sources. One such possibility presents itself if violation of the weak 
equivalence principle (VEP) occurs and the flavor eigenstates are not identical to the states 
that couple to gravity. The principle of equivalence is a cornerstone of Einstein's general 
theory of relativity. Normally such a premise would be considered sacrosanct, but so lit- 
tle has been experimentally tested for neutrinos that it may not be unreasonable to keep 
an open mind and check the validity of this principle for them. If this principle is in- 
deed violated then, as a consequence, the coupling of neutrinos to the gravitational field is 
nonuniversal. Under this circumstance, if the flavor eigenstates are linear superpositions of 



the gravitational eigenstates, VEP induced oscillations of neutrinos take place |12|. This 
does not require neutrinos to carry a non-zero mass. The important difference between 
this approach and the mass-mixing solution is manifested in the energy dependence of the 
survival probability. For a two-neutrino picture, the general expression for the survival 
probability for an initial Ve after propagation through a distance L is given by: 

Pee{E,, L) = l- sin2 29 sin2 {^\ , (1) 

where 9 is the rotation angle relating the gravitational eigenstate basis (i^i, 1^2) to the flavor 
basis {Uf.) i^xi X = //, r, or a sterile state). 

z/e = ui cos 9 + 1^2 sin 9] u^ = —vi sin 9 + 1^2 cos 9 (2) 

and A is the oscillation length, which for the VEP induced oscillation is: 



where E^ is the neutrino energy, (p the gravitational potential, and A/ = /i — /2 is a measure 
of the violation of equivalence principle, fi being the coupling strength of the gravitational 
eigenstates. In contrast, for mass- mixing vacuum oscillations A = {4:TrE^)/Am'^, where Am'^ 
is the mass square difference between two neutrino species. Thus, for the latter case \ (x E^ 
while for VEP A oc l/Ey. Due to the different energy dependences of the survival probability 
in the mass-mixing and the VEP alternatives, their predictions can be quite different. The 
phenomenological consequences of VEP-driven neutrino oscillations has attracted attention 



over the past decade [13| 



A completely unrelated situation which also leads to neutrino oscillations with A oc 1/Ei^ 
is a recently proposed picture of violation of special relativity (VSR) [jl4[. If special relativity 
is violated, the maximum attainable speed of a particle in vacuo need not universally be the 
speed of light c. In particular, if the maximum possible velocities of two types of neutrinos 
be vi and V2 and these velocity eigenstate neutrinos be related to the fe and i^x through 
a mixing angle 6 (see Eq.(|2|)) then the survival probability of a fg takes the same form as 
Eq.(||). In this case the expression for A is: 

where Av is the velocity difference for the neutrinos ui and 1/2- Comparing Eqs. (Q) and (|^ 
one finds that the energy dependence of the oscillation length is identical in the two case^ 
and the role of Av in the VSR case is the same as that of 4>Af in the VEP formalism. Here, 
we use the terminology of the VEP mechanism but the results can be taken over mutatis 
mutandis to the VSR situation. 

In this work we make a detailed examination of the VEP scenario in the light of the solar 
neutrino data. We consider the two possibilities of oscillation of the electron neutrino to (a) 
another active neutrino (t'^ or Vr), and (b) to a sterile (i.e., no weak interactions) neutrino. 
The Chlorine and Gallium experiments use radiochemical neutrino detection techniques 
and do not distinguish between oscillation of the z/g to an active neutrino or to a sterile 
one. The Kamiokande and SuperKamiokande experiments, on the other hand, are sensitive 
to other active neutrinos via the smaller neutral current contribution. Hence, the latter 
discriminates between the two alternatives of oscillation to active or sterile neutrinos. 

In section II, we present a brief summary of the ingredients that go into our analysis. 

In the next section we consider the total rates of solar neutrinos observed by the Chlorine, 
Gallex, SAGE, Kamiokande, and SuperKamiokande (1117-day data) experiments^ |^, ||, ^. 
We obtain the best-fit values of the VEP parameters and find that the goodness-of-fitfl 
(g.o.f ) is never high. We trace the origin of this to the difficulty within the VEP mechanism 
of simultaneously satisfying the Chlorine detector results and the very precise measurements 
of SK. We examine the effect of excluding either the Chlorine or the Kamiokande result in 



^It has been shown that inclusion of CPT-violating interactions in addition to Violation of Special Rela- 
tivity can lead to more general energy dependences involving 1/ E^, E^, and constant terms [|l5|] 
■^In the following, we refer to this data set as the 'total rates'. 
•^The goodness of fit gives the probability for the actual x^ to exceed Xn 
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the fits or of using the average of the two Gallium results rather than their individual mea- 
surements and find that the best-fit values are markedly different only in one situation, the 
exclusion of the Chlorine data from the analysis, when the quality of the fit is significantly 
improved. 

In section IV we turn to the 1117-day recoil electron energy spectrum from SK |I| and 
test the ability of the VEP model to account for the observations. In this case the quality 
of the fit is rather good. We examine a variant of the Standard Solar Model (SSM) in which 
the absolute normalization of the ®B neutrinos (Xb) is allowed to vary. In our analysis we 
have included the contribution from hep neutrinos and we have also examined the situation 
when the normalization of this flux (Xh) is different from the SSM value of unity. 

Section V deals with the simultaneous fitting of the total rates data and the SK 1117-day 
electron energy spectrum results. We find that good fits can be obtained, in this case, only 
if Xb and Xh are allowed to assume values different from the SSM stipulations. 

In section VI we turn to the expectations for SNO. We use the best fit values obtained 
from the above analyses to check how the predictions for the charged current deuteron 
disintegration and electron scattering are affected by VEP. 

We end in section VII with the conclusions. Similar work analysing solar neutrino data 



using the VEP formalism have been performed in the past |16, 17] and more recently in 



1 18, 19, 2C]. We compare our findings with these results. 



II Solar neutrinos and VEP 

The sun serves as a good neutrino factory. The fusion reactions that generate solar energy 
also produce neutrinos which are usually denoted by the different reactions {e.g., pp, Be, B, 
hep, etc.) from which they originate. The shape of the neutrino energy spectrum from any 
reaction is precisely known from weak interaction theory while the absolute normalization 
of the spectrum depends on the solar parameters like core temperature, opacity, etc. The 
pp neutrinos are the most copious but they are also of the lowest energy. Only the Gallium 
experiments (Gallex and SAGE) are sensitive to pp neutrinos and receive their dominant 
contributions from this process, in addition to smaller ones from the other reactions. At the 
other extreme, the Kamiokande and SuperKamiokande experiments are based on the water 
Cerenkov technique and can detect neutrinos of energy higher than about 5 MeV. Only the 
^B neutrinos (and to a small extent the hep neutrinos) contribute at these energies. The 
Chlorine experiment has a threshold of 0.8 MeV and mainly counts Be and ^B neutrinos. 
Thus the different detectors have probed different regions of the solar neutrino spectrum. 
Assuming that oscillations do occur during the passage of the neutrinos from their point 
of origin to the detector, the dependence of this phenomenon on the neutrino energy can 
be tested by the complementary information from the different experiments. From these 
results one should be in a position to check which of the vacuum, MSW, or VEP oscillations 
is preferred by the results. The SuperKamiokande collaboration has also presented the 
observed energy spectrum of the electron scattered by the incident neutrino. This spectrum 
is a direct probe of the energy dependence of neutrino oscillations. 



In this work we use the BP98 calculation [^] of the solar neutrino spectrum with INT 
normalization as the SSM reference. In addition to this, we explore the possibility of the 
absolute normalizations of the ^B- and /lep-neutrinos spectra being different from their BP98 
SSM predictions. If Xb and Xh denote the factors by which the absolute normalizations 
are multiplied, we use the data to find the best-fit values for these. In particular, we find 
that the fit to the SK recoil electron spectrum data is improved in a noteworthy manner 
when Xh is permitted to assume large values^. 

The weak equivalence principle requires the coupling of particles to the ambient gravita- 
tional potential <j) to be uniform, i.e., of the form f(j)E, where E is the particle energy, and / 
a universal coupling constant. If the latter varies from one neutrino species to another then 
that would constitute a violation of the equivalence principle. If /i 7^ /2 in a two-neutrino 
framework, then these states define a basis in the two-dimensional space which, in general, 
could be different from the flavor basis. The effect due to a small splitting A/ will manifest 
itself in the form of flavor oscillations, the wavelength going to infinity as A/ tends to zero. 
We follow the prevalent practice of choosing the gravitational potential, (p, to be a constant 
over the neutrino path. This is the case if the potential due to the Great Attractor Il2 
dominates over that due to the sun and other heavenly bodies in our neighborhood, 
such an event, writing AF = <pAf /2, the expression for the oscillation wavelength, Eq. 

becomes 

TT _ 6.20 X lO^i^m /lMeV\ 

" AF \ E, ) ■ 



A 



E^AF 
We choose the neutrinos to be massless. 



In 



(5) 



III Observed rates and VEP oscillation of neutrinos 

The data used for the x^ analysis of total rates are given in Table 1. They are from 
the Chlorine experiment at Homestake, the two Gallium experiments, Gallex and SAGE, 
the water Cerenkov detector Kamiokande, [^, and the latest 1117-day data from Super- 
Kamiokande M. 



Experiment 


Chlorine 


Gallium 


Kamiokande 


Super 
Kamiokande 


Gallex 


SAGE 


Observed Kate 
BP98 Prediction 


0.33 ± 0.029 


0.60 ±0.06 


0.52 ± 0.06 


0.54 ±0.07 


0.465 ±0.015 



Table 1: The ratio of the observed solar neutrino rates to the corresponding BP98 SSM predictions 
used in this analysis. The results are from Refs. y| and |[|]. 



The definition of x^ used for this analysis is: 



X 



J2 (i^*/^ _ Fr) {a-;") {Fj'^ - Ff ^) . 



(6) 



«J=1,5 



''See, however, remarks in the penultimate paragraph of section VII on the fits to the separate day and 
night SK electron spectra measurements. 



Here F^ = j"^ Jj'pP^^ where ^ is th (for the theoretical prediction) or exp (for the exper- 
imental value) and Ti is the total rate in the z-th experiment. F^^^ is taken from Table 
1. The error matrix aij contains the experimental errors, the theoretical errors and their 
correlations. The theoretical errors have contributions which originate from uncertainties 
in the detector cross-sections as well as from astrophysics |23]. The off-diagonal elements 
in the error matrix come through the latter. In estimating the astrophysical contribution, 
the uncertainties in the spectrum due to the input parameters are taken from pT| . 

In the presence of neutrino conversions, the detection rate on earth for the radiochemical 
Chlorine and Gallium experiments is predicted to be: 



-,th 



E 



Eth 



Xc,<l)a{Eu)Oi{Ey) < Pee{Ey) >« dEy, 



(7) 



where (Ti{Ey) is the neutrino capture cross-section for the i-th detector |24] and Et^ the 
neutrino threshold energy for detection. (j)a{Ey) stands for the neutrino spectrum for the 
Q-th source [p4| and Xa is an overall normalization factor for this spectrum such that X^ = 1 
corresponds to the SSM. The sum is over all the individual neutrino sources. < Pee{Eu) >« 
is the neutrino survival probability for the a-th source averaged over the distribution of 
neutrino production regions in the sun, 



< PeeiEr.) >a= J dr Pee{E^,r)^a{r). 



(8) 



<I>Q,(r) is a normalized function which gives the probability of the a-th reaction occuring at 
a distance r from the center of the sun and Pee{E^, r) is obtained by averaging the survival 
probability over a year taking the eccenticrity of the earth's orbit into account, i.e., 



P,e{E,,r) = ^ f dt 
1 Jo 

where R{t) is the sun-earth distance given by. 

R{t) = Ro 



^ r o ^ o (TTR(t) 

1 - W 29 sin^ ' ^ ^ 



R{t) 



l-ecos(27r-) 



(9) 



(10) 



Here, Rq = 1.49 x 10 m is the mean Sun-Earth distance and e = 0.0167 is the ellipticity 
of the earth's orbit, t is the time of the year at which the solar neutrino flux is measured 
and T is 1 year. 

The theoretical prediction according to the BP98 Standard Solar Model, Tj , is 
obtained by setting the survival probability as 1 in the above. 

For the water Cerenkov detectors Kamiokande and SuperKamiokande, in the case of 
oscillation to active neutrinos one has to take into account the contributions to the signal 
from the produced f^ or z^,- via the neutral current interactions. 
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dE^ / dEx 
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< PeeiE.) >a ^+ < Pe,{E.) > '^'"'^ 



dEr 



dEr 



(11) 



The second term in the bracket is absent if oscillation to sterile neutrinos is under consider- 
ation. Et and Ea denote the true and apparent electron energies respectively. Ej-^.^ and 
^Tmax ^^^ determined by kinematics. p{Ea,Et) is the energy resolution function for which 
we use the expression given in |^. Ea^,^ is 7.5 (5.5) MeV for the calculation of the total 
rate at Kamiokande (SuperKamiokande). The differential cross-section for the production 
of an electron with true relativistic energy Et, j;§-, is obtained from standard electroweak 
theory. 

Now we discuss the results obtained from a x^ minimization analysis of the data within 
the VEP picture. For these fits we have set X^ = 1 for all a; i.e., we take the normalizations 
of the solar spectra at their SSM values. The best-fit parameters, Xmin/i^^S^^^ of freedom), 
and the goodness of fit values are presented in Table 2. Both possibilities of oscillation of 
the I'e to an active or a sterile neutrino have been considered. To gauge the impact of the 
different experiments on the best-fit values of the VEP parameters, we have first fitted all 
the five data given in Table 1. We have then repeated the procedure excluding the Chlorine 
or the Kamiokande results. We also examine the effect of using the average of the result of 
the two Gallium measurements (0.57ib 0.054) rather than the two individual ones. It is seen 
from Table 2 that in the different alternatives the best-fit values of the parameters are all 
rather close excepting for the case where the Chlorine result is excluded from the analysis. 
The fit improves significantly in the latter case; but even here the goodness of fit for active 
(sterile) neutrinos is still just 36% (58.5%). It is seen that maximal mixing is preferred and 
the VEP oscillation length is (see Eq. (^) ~ 3.4 x 10^^ m, comparable with the earth-sun 
distance. This is reminiscent of the mass-mixing vacuum oscillation solution to the solar 
neutrino problem. 



Neutrino 
Type 


Set 


Fitted 
Experiments 


sin^ 29 


AF 

(10-24) 


xLn/d-o.f 


g.o.f 

(%) 


Active 


la 


CI, Gallex, SAGE, K, SK 


1.0 


1.80 


4.54/3 


20.84 


lb 


Gallex, SAGE, K, SK 


0.85 


4.57 


2.04/2 


36.05 


Ic 


CI, Gallex, SAGE, SK 


1.0 


1.80 


3.50/2 


17.41 


Id 


CI, (Ga)av, K, SK 


1.0 


1.80 


3.72/2 


15.59 


le 


CI, (Ga)av, SK 


1.0 


1.80 


2.67/1 


10.20 


Sterile 


2a 


CI, Gallex, SAGE, K, SK 


1.0 


1.84 


5.89/3 


11.73 


2b 


Gallex, SAGE, K, SK 


0.84 


4.53 


1.07/2 


58.51 


2c 


CI, Gallex, SAGE, SK 


1.0 


1.84 


4.58/2 


10.12 


2d 


CI, (Ga)av, K, SK 


1.0 


1.84 


5.08/2 


7.91 


2e 


CI, (Ga)av, SK 


1.0 


1.84 


3.77/1 


5.23 



Table 2: The best-fit values of the parameters, sin^ 29, AF, x 
rates of the different experiments. 



2 

min' 



and the g.o.f. for fits to the total 



As is seen from Table 2, the goodness of fits are rather poor for all the above cases. 
In order to trace the origin of this result, we present in Table 3 the total rates for the 
different experiments obtained using the best-fit values of the VEP parameters presented 



in Table 2. These should be compared with the experimental data in Table 1. Several 
points are noteworthy. The rates for the different experiments are not sensitive to the 
small changes in the input VEP parameters; the numbers for the Gallium experiments and 
Kamiokande are always within la of the experimental value but for Chlorine the deviation 
is 3a. The fit to SuperKamiokande is always bad {8a) irrespective of which experiments are 
excluded from the analysis excepting for the singular case where the Chlorine rate is left 
out when a 2a fit is obtained. This is a reflection of the very precise nature of the present 
SK data and the inability of the VEP mechanism to simultaneously reproduce the varying 
degrees of suppression seen in experiments with different energy thresholds. In particular, 
the Chlorine result is seen to be especially problematic in this respect. It can be surmised 
from the Chlorine-excluded analyses - Table 3 (sets lb, 2b) - that a simultaneous good fit 
to all data would be obtained if the suppression in the Chlorine experiment had been less 
than that in the Gallium experiments (lower threshold) as well as in the water Cerenkov 
ones (higher threshold) , indicating that for these fit-values of the VEP parameters the "^Be 
neutrinos are less suppressed than the other solar neutrinos. 



Neutrino 
Type 


Set 


Fitted 
Experiments 


CI 


Ga 


K 


SK 


Active 


la 


CI, Gallex, SAGE, K, SK 


0.436 


0.572 


0.600 


0.600 


lb 


Gallex, SAGE, K, SK 


0.655* 


0.549 


0.508 


0.497 


Ic 


CI, Gallex, SAGE, SK 


0.436 


0.572 


0.600* 


0.600 


Id 


CI, (Ga)av, K, SK 


0.436 


0.574 


0.600 


0.600 


le 


CI, (Ga)av, SK 


0.436 


0.574 


0.600* 


0.600 


Sterile 


2a 


CI, Gallex, SAGE, K, SK 


0.435 


0.565 


0.525 


0.525 


2b 


Gallex, SAGE, K, SK 


0.657* 


0.555 


0.512 


0.468 


2c 


CI, Gallex, SAGE, SK 


0.435 


0.565 


0.525* 


0.525 


2d 


CI, (Ga)av, K, SK 


0.435 


0.565 


0.525 


0.525 


2e 


CI, (Ga)av, SK 


0.435 


0.565 


0.525* 


0.525 



Table 3: The total rate predictions for the different experiments obtained by using the best-fit values 
of the VEP parameters presented in Table 2. Predictions for experiments not included in the fit 
are marked with an asterisk. 



We show in Fig. 1, the allowed region for the parameters sin^ 20 and AF at 90% 
C.L. for the (la) active and (lb) sterile cases obtained by fitting the total rates of all 
five experiments, i.e., Chlorine, Gallex, SAGE, Kamiokande, and SuperKamiokande. The 
confidence level has been fixed with respect to the global minimum. The best fit points 
have been indicated. It is seen from Fig. 1 that the nature of the allowed regions for the 



active and sterile cases are roughly similar, an observation made by Gago et al. [19|. The 
sterile alternative is more restrictive. 



IV SuperKamiokande Recoil electron energy spectrum and 
VEP 



We now turn to a fit to the recoil electron energy spectrum as seen at SuperKamiokande. 
We use the 1117-day data for this analysis. The SK results have been presented in the form 
of number of events in 17 electron recoil energy bins of width 0.5 MeV in the range 5.5 MeV 
to 14 MeV and an 18th bin which covers the events in the range 14 to 20 MeV [|l|, p6[ . 

With the SuperKamiokande threshold in the 5 MeV range, only the ^B and hep solar 
neutrinos contribute to the signal. In addition to the SSM, for this analysis we have ex- 
amined models in which the normalizations of both the ^B and hep neutrino flux (Xb and 
Xh, respectively - normalized to unity for the SSM) are allowed to vary arbitrarily. In this 
case, x^ is defined as 



X 



i,j=l,18 



(i?f 



Rr") 



a. 



ij )sp 



(« 



th 
J 



RfP 



(12) 



Ri = Sf /SP^^^ with ^ being th or exp as before and Si standing for the number of events 
in the i-th energy bin. The theoretical prediction is given by eq. (^) but the integration 
over the apparent (i.e., measured) energy will now be over each bin. The error matrix aij 
used by us is Q 



{-' 



ijjsp 



Sij [ cr, 



2 ,2 

i,stat "T ^' 



i,uncorr. 



~r CTiexp^jjCxp ~r ^i^cal^j,cah 



(13) 



where we have included the statistical error, the uncorrelated systematic errors and the 



energy-bin-correlated experimental errors |27] as well as those from the calculation of the 
shape of the expected spectrum [^]. Since we allow the normalizations of the ®B and hep 
fluxes to vary, we do not include their astrophysical uncertainties separately. The results 
are presented in Table 4. 



Neutrino 
Type 


Set 


sin^ 26 


AF 

(10-24) 


Xb 


Xh 


Xmin,/d-0.f. 


g.o.f 

(%) 


Active 


3a 


1.0 


1.97 X 10^ 


0.79 


15.07 


8.83/14 


84.2 


3b 


0.38 


0.22 


0.65 


1.0 

(fixed) 


10.04/15 


81.7 


3c 


0.69 


0.23 


1.0 

(fixed) 


-2.20 


11.59/15 


70.97 


3d 


0.68 


0.23 


1.0 

(fixed) 


1.0 

(fixed) 


11.66/16 


76.7 


Sterile 


4a 


1.0 


39.33 


1.03 


11.52 


8.78/14 


84.48 


4b 


0.34 


0.22 


0.67 


1.0 

(fixed) 


10.01/15 


81.9 


4c 


0.98 


39.33 


1.0 

(fixed) 


11.94 


8.81/15 


88.71 


4d 


0.58 


0.23 


1.0 

(fixed) 


1.0 

(fixed) 


11.34/16 


78.8 



Table 4: The best-fit values of the parameters, sin 29, AF, Xb, X^, Xmiw ^'^'^ "the g.o.f. for fits to 
the scattered electron spectrum at SK. 

The large values of X^ obtained from the fits with the best g.o.f. (3a and 4b) indicate 
that an increased number of high energy hep neutrinos yields a better fit to the highest 
energy bin of the spectrum observed at SK. It may be noteworthy that a rise in the observed 
electron energy spectrum at the high energy end seen in the earlier 825-day data has become 
less prominent in the latest 1117-day sample. A further softening will bring Xh closer to 
the SSM prediction^]. 



Neutrino 
Type 


Set 


sin^ 2d 


AF 

(10-24) 


Xb 


Xh 


CI 


Ga 


K 


SK 


Active 


3a 


1.0 


1.97 
xl02 


0.79 


15.07 


0.412 


0.463 


0.471 


0.466 


3b 


0.38 


0.22 


0.65 


1.0 

(fixed) 


0.585 


0.947 


0.481 


0.488 


3c 


0.69 


0.23 


1.0 
(fixed) 


-2.20 


0.560 


0.944 


0.497 


0.503 


3d 


0.68 


0.23 


1.0 

(fixed) 


1.0 

(fixed) 


0.562 


0.944 


0.507 


0.513 


Sterile 


4a 


1.0 


39.33 


1.03 


11.52 


0.562 


0.538 


0.479 


0.477 


4b 


0.34 


0.22 


0.67 


1.0 

(fixed) 


0.612 


0.951 


0.485 


0.492 


4c 


0.98 


39.33 


1.0 
(fixed) 


11.94 


0.563 


0.548 


0.477 


0.476 


4d 


0.58 


0.23 


1.0 

(fixed) 


1.0 

(fixed) 


0.630 


0.953 


0.507 


0.509 



Table 5: The calculated values of the rates for the different experiments using the best-fit values of the 
parameters, sin'^ 26*, AF, Xb, and X^, from fits to the scattered electron spectrum at SK. 

It may be of interest to check how the best-fit values of the VEP parameters fare when 
confronted with the total rates data. For this purpose, we use the four sets of best-fit values 
of parameters from Table 4 and use them to compute the total rates from the different 
experiments. These results are shown in Table 5. The latter may be compared to the values 
obtained using fits to the total rates themselves (see Table 3). Using Table 1, it is seen 
that the predicted rate for SK is in rather close agreement with the measured value, as is to 
be expected, but for the Chlorine and Gallium experiments the deviations are very large. 
Notice, in particular, that in several cases only a very tiny suppression is obtained for the 
Gallium experiments, though for the fits to the spectrum with the best g.o.f. - fits 3a and 
4c (see Table 4) ~ this is not the case. This underscores a basic incompatibility, within the 
VEP picture, of the results from the different experiments. 

^See the remarks on fits to the SK separate day and night spectra measurements in section VII 
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In Fig. 2 we show the 90% C.L. allowed region in the sin 26-AF plane obtained by 
fitting the scattered electron spectrum for the Xb = X^ = 1 case (area enclosed by the 
solid lines). The best-fit point is marked with a dark dot. Also shown is the 90% C.L. 
disallowed region when Xb is permitted to float arbitrarily but keeping X^ fixed at unity 
(area enclosed by the broken lines). The best-fit point in this case has been indicated by a 
X sign. Both (a) active and (b) sterile neutrino alternatives have been displayed. 



V Combined analysis of rates and spectrum in the VEP sce- 
nario 

In the previous sections we have examined the solar neutrino total rates and the SK scat- 
tered electron energy spectrum within the VEP oscillation framework. For the latter, we 
found a good fit using SSM input parameters and even better ones when the absolute nor- 
malizations of the *B and/or hep fiuxes are allowed to vary (see Table 4). The rates fit, 
on the other hand, was not satisfactory and improved somewhat if the results from the 
Chlorine experiment was left out from the x^ analysis. In this section we make a combined 
fit to the rates and spectrum data. 

For the combined analysis of rate and spectrum, we take the rate measurements from 
the Chlorine experiment, the average of the Gallex and SAGE results, and the SK 1117-day 
rate data. To obtain the x^ corresponding to any value of the input parameters, separate 
values of x^ calculated for the total rates and the spectrum data are added and then the 
total is minimized. The best-fit values for VEP oscillations to active and sterile neutrinos 
are given in Table 6. The number of degrees of freedom for this case is (18-1-3 - 4)=17. 
Notice that good fits (cases 5a and 6a) can be obtained only when both Xb and X^ are 
allowed to assume values different from their SSM requirements. 



Neutrino 
Type 


Set 


sin^ 26 


AF 

(10-24) 


Xb 


Xh 


xLn/d-o.f. 


g.o.f 

(%) 


Active 


5a 


1.0 


1.72 


0.79 


23.69 


11.51/17 


82.87 


5b 


1.0 


1.68 


0.79 


1.0 

(fixed) 


12.96/18 


79.38 


5c 


1.0 


19.26 


1.0 

(fixed) 


33.57 


50.79/15 


8.9 xlO-4 


5d 


1.0 


18.75 


1.0 

(fixed) 


1.0 

(fixed) 


70.44/16 


8.37 xlO-'' 


Sterile 


6a 


1.0 


1.77 


0.91 


34.55 


14.04/17 


66.45 


6b 


1.0 


1.86 


0.91 


1.0 

(fixed) 


16.28/18 


57.28 


6c 


1.0 


1.86 


1.0 

(fixed) 


36.87 


26.90/15 


2.95 


6d 


1.0 


1.89 


1.0 

(fixed) 


1.0 

(fixed) 


29.84/16 


1.88 
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Table 6: The best-fit values of the parameters, sin 29, AF, Xb, Xh, Xmiw ^'^'^ "the g.o.f. for fits to 
the total rates as well as the scattered electron spectrum at SK. 

In exhibiting the 90% C.L. allowed regions in parameter space from this combined fit, 
we present only the results for the Xh = 1,Xb arbitrary case. As already noted, if the 
SSM normalization for the ^B flux is used {Xb = 1) then the goodness-of-fits are poor. 
These cases are not pursued further. In Fig. 3 we present the 90% C.L. allowed regions in 
the sin 29 - AF plane from a combined analysis of total rates and spectrum. The best-fit 
points are marked by a x sign. Both (a) active and (b) sterile neutrino alternatives have 
been considered. It is noteworthy that the allowed regions have a considerable overlap with 
that in Fig. 1 obtained from a fit to the total rates alone. 

VI Predictions for SNO 

The analyses described in the previous sections yield sets of best-fit values for the parameters 
sin^20, AF, Xb and X^- In this section we examine what these best-fit values imply for 
the Sudbury Neutrino Observatory (SNO) experiment ||29t| . 

At SNO the neutrinos are detected by three processes, namely, (a) charge current (CC) 
break up of the deuteron, (b) electron scattering by the neutrino, and (c) neutral current 
(NC) break up of the deuteron. 

(14) 
(15) 
(16) 

For the scattering ([l^ ) and NC ( |l6|) reactions, v stands for any active neutrino; it has to 
be borne in mind that for the fi and r fiavours only the Z exchange contribution is present 
for the former reaction while for the I'e there is an additional (dominant) piece from W- 
exchange. For the CC reaction and for scattering, the electrons are detected by the emitted 
Cerenkov radiation and hence their energy spectrum is directly measured as at SK. For the 
NC reaction, on the other hand, only calorimetric measurment is possible. At present, data 
is being taken for the first two processes. 

For the scattered electrons, the formalism for theoretical prediction is similar to that 



V'e + d - 


->■ p + p + e~ 


(CC reaction) 


f + e~ - 


-> V + e~ 


(scattering) , 


V + d - 


->■ V + p + n 


(NC reaction). 



for SK (see eq. (11)) except for the fact that here in SNO, the detector fluid is IK ton of 
D2O instead of 32K ton of water in SK. 

For the CC interactions, the theoretical predictions for the electron energy spectrum at 
SNO is given by 

Tsho = E / dE, dEr / dEA p{Ea, Et) X^ MEu) < Pee{E,) >„ a^^ . 

(17) 
The charge current scattering cross-section at SNO, a , is found from Ref. |3^. In our 
calculations, we use the resolution function p{Ea,Et) for SNO as given in ||25|. We have 
set Ea^^ to 5.0 MeV. 
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In Fig. 4 we show the signal expected at SNO due to (a) v — e scattering and (b) 
the CC reaction if the VEP oscillation is operative. We choose two typical cases for the 
VEP parameters: (i) the best-fit values obtained from the fit to the spectrum using the 
SSM inputs {Xb = Xh = 1), and (ii) the best-fit parameters from the combined rate and 
spectrum fit when both Xb and X^ are allowed to vary. The SSM expectation in the 
absence of VEP transitions is also shown. 

Figs 5(a) and 5(b) are similar to Figs. 4(a) and 4(b) excepting that we present the expec- 
tations at SNO normalized by the SSM predictions. This form may be of more convenience 
for comparison with the experimental results. 

SNO will subsequently also detect neutrinos via the NC reaction (16). Since all active 



neutrinos register in the NC reaction with the same strength, this reaction provides a means 
to distinguish oscillations to active neutrinos from those to sterile ones. The ratio of the 
NC and CC rates, Rnc and Rcc^ is somewhat less sensitive to theoretical uncertainties 
than the rates themselves. Therefore, for the purpose of illustration, we present in Table 7 
the values for the ratios 

RsNO = 15 , (18) 

Rcc 

Q n /TjBP98 Rnc/RnC 
^SNO - RSNO/RsNO - -^ /RBP98 ' 



For the NC cross-sections we use Ref. |30]. Results are presented for the best-fit values of 
the parameters obtained from the total rates, the SK scattered electron spectrum, and the 
combined rates and spectrum data presented in Tables 2, 4, and 6. 

It is seen from Table 7 that the variables Rsno and Ssno are particularly effective in 
distinguishing VEP oscillations to active neutrinos from those to a sterile neutrino. For ac- 
tive neutrinos Rsno {Ssno) varies between 0.44 - 0.74 (1.15 - 1.92) while the corresponding 
range for sterile neutrinos is 0.29 - 0.30 (0.77 - 0.79). 
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Neutrino 
Type 


Set 


sin^ 26 


AF 

(10-24) 


Xb 


Xh 


RSNO 


SsNO 


Active 


la,c,d,e 


1.0 


1.80 


1.0 


1.00 


0.610 


1.57 


lb 


0.85 


4.57 


1.0 


1.00 


0.740 


1.92 


3a 


1.0 


1.97 X 10^ 


0.79 


15.07 


0.581 


1.50 


3b 


0.38 


0.22 


0.65 


1.0 

(fixed) 


0.443 


1.15 


3c 


0.69 


0.23 


1.0 

(fixed) 


-2.20 


0.724 


1.87 


3d 


0.68 


0.23 


1.0 

(fixed) 


1.0 

(fixed) 


0.712 


1.89 


5a 


1.0 


1.72 


0.79 


23.69 


0.583 


1.51 


Sterile 


2a,c,d,e 


1.0 


1.84 


1.0 


1.00 


0.305 


0.789 


2b 


0.84 


4.53 


1.0 


1.00 


0.303 


0.784 


4a 


1.0 


39.33 


1.03 


11.52 


0.300 


0.776 


4b 


0.34 


0.22 


0.67 


1.0 

(fixed) 


0.299 


0.774 


4c 


0.98 


39.33 


1.0 

(fixed) 


11.94 


0.299 


0.775 


4d 


0.58 


0.23 


1.0 

(fixed) 


1.0 

(fixed) 


0.302 


0.781 


6a 


1.0 


1.77 


0.91 


34.55 


0.302 


0.781 



Table 7: The ratio of the NC rate to the CC rate at SNO, Rsno, and the same ratio normalized to the 
BP98 SSM, SsNO for the best-fit values of the parameters, sin^ 29, AF, Xb, X^ obtained from 
fits to the total rates and the SK scattered electron spectrum. The sets 1 and 2 are from Table 2, 
3 and 4 from Table 4, and 5 and 6 from Table 6. 

VII Discussions and Conclusions 

In this work we have made a detailed examination of the viability of the VEP oscillation 
mechanism in the light of the solar neutrino data. The parameters in the VEP formalism are 
AF, a measure of the violation of the weak equivalence principle, and 9, the mixing angle 
relating the flavor basis of neutrinos to the gravitational basis. The data we have included 
in the analysis come from the radiochemical Chlorine experiment, the similar Gallium- 
based Gallex and SAGE collaborations, and the Cerenkov technique reliant Kamiokande 
and SuperKamiokande (1117-day data) facilities. We have found that good flts (x^/d.o.f. 
= 11.66/16 for active neutrinos) can be obtained to the scattered electron energy spectrum 
seen at SK within the VEP mechanism using SSM inputs. The fits are even better in models 
in which the absolute normalizations of the ^B and/or hep fluxes are allowed to vary. The 
preferred oscillation wavelength is comparable to the earth-sun distance and is reminiscent 
of the mass-mixing vacuum oscillation solution. In contrast to the flt to the spectral data, 
the fit to the observed rates is poor but improves somewhat (x^/d.o.f. = 2.04/2 for active 
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neutrinos) if the Chlorine experiment is left out of the fit. A combined fit to the rates and 
spectrum has also been performed and it is found that within the SSM the goodness is 
not satisfactory. However, in a model in which the absolute normalizations of the '^B and 
hep fluxes are allowed to vary, the fit is much improved (x^/d.o.f. = 11.51/17 for active 
neutrinos with Xb = 0.79 and Xh = 23.69). The best-fit values of AF (= 1.7 xlQ-^^) and 
sin^ 20 (= 1) are closer to those obtained from the fits to the rates alone. We find that the 
90% C.L. allowed regions are broadly similar whether VEP transitions occur to active or 
sterile neutrinos. 

It is of interest to compare the results that we have obtained with those of other recent 
analyses of the solar neutrino problem in the VEP picture |19, 20 1. For the fit to the total 



rates, our results agree almost exactly with those of [20| and are within 5% of those of |19|. 
The earlier analyses used the 825-day SK results and we must therefore conclude that the 
newer data do not make a significant impact on this fit. In contrast, for the fits to the SK 
recoil electron energy spectrum, our best-fit values of /S.F and sin^ 26 are both somewhat 
smaller than those of |jl^, |2y] . We have checked that if we use the 825-day SK data there is 
better agreement. We attribute the difference to the reduction of the higher energy events 
in the newer data. For the combined fit to the rate and spectrum, we are in good agreement 
with [|l^; the goodness-of-fit has improved with the newer data. 

During the passage of the neutrinos from their point of production to the solar surface, 
interactions with the ambient matter, responsible for the MSW effect, become important. 
Apart from a neutral current contribution which affects the masses of all active neutrino 
species identically (and is therefore irrelevant for this discussion), there is a contribution to 
the electron neutrino mass niMSW — V^Gpn^ir), Gp being the Fermi coupling and nf.{r) 
the number density of electrons at a distance r from the center of the sun, due to charge 
current interactions. In the presence of VEP and this MSW contribution, the effective 
neutrino mass matrix in flavor space takes the form 

E^I^F cos2e -2^/2GFne{r) E„AFsm2e 
E^AFsm29E^ -E^AFcos2e 



2 



(19) 



where we have dropped an irrelevant part proportional to the identity matrix. The MSW 
contribution in ( p!9| ) inside the sun turns out to be several orders of magnitude larger than 
the terms due to VEP that we have discussed in this work. Recall that we have assumed 
the neutrinos to be massless. For the maximal mixing case (i.e., sin 20 = 1), which we 
have found for the best-fit solutions, there is no resonance effect and, in fact, till such time 
that the neutrino emerges from the sun, the MSW contribution controls the masses in (|l9|). 
Inside the sun, the Ve is, therefore, a mass eigenstate to a very good approximation. The 
effect of VEP oscillations begins to manifest itself only from then onwards. For oscillation 
to a sterile neutrino, the neutral current contribution to the Ve mass also becomes relevant. 
It is of the same order as the charged current piece and the resultant effect in this case is 
much the same as that for the active neutrino alternative. There is a different region in the 
sin 29, AF parameter space (larger AF, non- maximal mixing) where other solutions have 
been found with the MSW effect playing an important role |l^, [l^p|. 

''It has, however, been recently demonstrated that these solutions are inconsistent with the results from 
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SuperKamiokande has now obtained the recoil electron spectra separately for day and 



night runs using their 1117-day data [31|. The oscillation wavelengths corresponding to the 
best-fit VEP parameters that we have found are comparable to the earth-sun distance (1.49 
xlO^^m). The extra distance travelled by neutrinos through the earth during the night 
runs is insignificant compared to this. This distance can be readily estimated. Using the 
latitude of SuperKamiokande (36.4°) and the obliquity of the ecliptic (23.5°), we estimate 
the maximum extra distance travelled during the night run to be 6.40 xlO^m. To confirm 
the expectation that fits to the separate day and night spectra give results which are not 
much different, we performed a ^^-minimization on this data for oscillation to an active 
neutrino. When Xb and X^ were held fixed at their SSM values of unity, the best fit 
corresponds to sin 29 = 0.71 and AF = 0.21 x 10"^"^. For comparison, the corresponding 
values from Table 4, set 3d are sin^ 261 = 0.68 and AF = 0.23 x lO'^^. If Xb and Xh 
are allowed to vary then we get sin^ 2(9 = 1.0 (1.0), AF = 1.97 x lO'^^ (1.97 x 10"^^), 
Xb = 0.78 (0.79), and Xh = 0.74 x 10"^ (15.07). In the parantheses we have given the 
values presented earlier in Table 4, set 3a. Notice that there is no significant change barring 
the very small value of Xh- The Xh dependence of x^ is slow and this indicates that the 
day-night scattered electron spectra is consistent with no hep neutrinos. 

The VEP oscillation wavelengths favored by the data are comparable to the distance 
of the sun to the earth. Therefore, seasonal effects are to be expected if this mechanism 
is responsible for the solar neutrino deficit. However, the data from SK is still not of very 
high statistics and we have left this analysis to a future work. 
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Figure Captions 

Fig. 1. The 90% C.L. allowed region in the sin^ 29 - AF plane from an analysis of the total 
rates seen by the Chlorine, Gallex, SAGE, Kamiokande, and SuperKamiokande (1117-day 
data) detectors assuming conversion due to VEP to (a) active neutrinos and (b) sterile 
neutrinos. The best fit points have been indicated. 

Fig. 2. The region bounded by the solid lines is the 90% C.L. allowed region in the sin^ 26 
- AF plane from the fitting of the 1117-day SK recoil electron spectrum data with the 
SSM flux normalizations (Xb = X^ = 1). The best fit point is marked with a dark dot. 
The best-fit point for the case Xh = 1,Xb arbitrary, is indicated by a 'x' sign and the 
corresponding 90% C.L. disallowed region is enclosed by the broken lines. The two panels 
correspond to (a) active neutrinos and (b) sterile neutrinos. 

Fig. 3. The 90% C.L. allowed region in the sin^ 29 - AF plane from the combined analysis 
of total rates and 1117-day SK recoil electron spectrum data for conversion due to VEP. 
The hep flux normalization is held fixed at the SSM value (X^ = 1) but the ^B flux nor- 
malization {Xb) is allowed to vary. The best fit points are indicated. Panel (a) represents 
the case of conversion to an active neutrino while panel (b) corresponds to conversion to a 
sterile neutrino. 

Fig. 4. The number of events expected per day at SNO due to (a) i' — e scattering, and (b) 
the charged current U(, — d reaction. Here E^. stands for the observed electron energy. Re- 
sults are shown for the best-fit values of the parameters obtained from fits to the spectrum 
and jointly to the rates and the spectrum (see text for more details). In both panels the 
solid (large dashed) line corresponds to the best-fit parameters obtained from a fit to the 
electron energy spectrum for VEP oscillation to an active (sterile) neutrino. The dot-dashed 
(dotted) lines similarly correspond to the best combined fits to the rates and spectrum for 
VEP oscillations to an active (sterile) neutrino. The BP98 SSM expectations are the small 
dashed curves. 

Fig. 5. The ratio of the predicted number of events from VEP to the SSM expectation as 
a function of the electron energy for SNO due to (a) v — e scattering, and (b) the charged 
current V(, — d reaction (see text for more details). Here E^ stands for the observed electron 
energy. In both panels the solid (large dashed) line corresponds to the best-fit parameters 
obtained from a fit to the electron energy spectrum for VEP oscillation to an active (sterile) 
neutrino. The dot-dashed (dotted) lines similarly correspond to the best combined fits to 
the rates and spectrum for VEP oscillations to an active (sterile) neutrino. 
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